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ABSTRACT: The ground state density matrix for a massless free field is 
traced over the degrees of freedom residing inside an imaginary sphere; the resulting 
entropy is shown to be proportional to the area (and not the volume) of the sphere. 
Possible connections with the physics of black holes are discussed. 
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A free, massless, scalar, quantum field (which could just as well represent, say, the 
acoustic modes of a crystal, or any other three-dimensional system with dispersion rela- 
tion 00 = c\k \ ) is in its nondegenerate ground state, |0). We form the ground state density 
matrix, pq = |0)(0|, and trace over the field degrees of freedom located inside an imaginary 
sphere of radius R. The resulting density matrix, p ut 5 depends only on the degrees of free- 
dom outside the sphere. We now compute the associated entropy, S = — Tr p ut log Pout- 
How does S depend oni?? 

Entropy is usually an extensive quantity, so we might expect that S ~ R?. However, 
this is not likely to be correct, as can be seen from the following argument. Consider tracing 
over the outside degrees of freedom instead, to produce a density matrix pj n which depends 
only on the inside degrees of freedom. If we now compute S = — Tr p- m log p- m , we would 
expect that S scales like the volume outside the sphere. However, it is straightforward to 
show that p- m and p ou t have the same eigenvalues (with extra zeroes for the larger, if they 
have different rank), so that in fact S = S [1]. This indicates that S should depend only 
on properties which are shared by the two regions (inside and outside the sphere). The 
one feature they have in common is their shared boundary, so it is reasonable to expect 
that S depends only on the area of this boundary, A = AuR?. S is dimensionless, so to get 
a nontrivial dependence of S on A requires another dimensionful parameter. We have two 
at hand: the ultraviolet cutoff M and the infrared cutoff p, both of which are necessary 
to give a precise definition of the theory. (For a crystal, M would be the inverse atomic 
spacing, and p the inverse linear size, in units with h = c = 1.) Physics in the interior 
region should be independent of p, which indicates that perhaps S will be as well. We 
therefore expect that 5 is some function of M 2 A. 

In fact, as will be shown below, S = kM^A, where k is a numerical constant which 
depends only on the precise definition of M that we adopt. 

This result bears a striking similarity to the formula for the intrinsic entropy of a black 
hole, >Sbh = \M^ x A, where M P1 is the Planck mass and A is the surface area of the horizon 
of the black hole [2]. The links in the chain of reasoning establishing this formula are 
remarkably diverse, involving, in turn, classical geometry, thermodynamic analogies, and 
quantum field theory in curved space. The result is thus rather mysterious. In particular, 
we would like to know whether or not S B n has anything to do with the number of quantum 
states accessible to the black hole. 

As a black hole evaporates and shrinks, it produces Hawking radiation whose entropy, 
S'hr, can be computed by standard methods of statistical mechanics. One finds, after the 
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black hole has shrunk to negligible size, that S'hr is a number of order one (depending on 
the masses and spins of the elementary particles) times the original black hole entropy [3] . 
This calculation of S HR is done by counting quantum states, and the fact that S Bii ~ S HR 
lends support to the idea that S BH should also be related to a counting of quantum states. 
It is then tempting think of the horizon as a kind of membrane [4], with approximately 
one degree of freedom per Planck area. However, in classical general relativity, the horizon 
does not appear to be a special place to a nearby observer, so it is hard to see why it 
should behave as an object with local dynamics. The new result quoted above indicates 
that S ~ A is a much more general formula than has heretofore been realized. It shows 
that the amount of missing information represented by S Bli is about right, in the sense 
that we would get the same answer in the vacuum of flat space if we did not permit 
ourselves access to the interior of a sphere with surface area A, and set the ultraviolet 
cutoff to be of order M P1 (perfectly reasonable for comparison with a quantum theory that 
includes gravity). Furthermore, getting S ~ A clearly does not require the boundary of 
the inaccessible region to be dynamical, since in our case it is entirely imaginary. 

To establish that S = kM 2 A for the problem at hand, let us begin with the simplest 
possible version of it: two coupled harmonic oscillators, with hamiltonian 

H = \[pI +P2 + k (4 + 4) + h(xi -x 2 f] ■ (1) 

The normalized ground state wave function is 

ip (x 1 ,x 2 ) = 7r~ 1/2 (u+U-) 1/4: exp[-(u; + x\+uj-X 2 _)/2] , (2) 

where x± = (x\ ± x 2 )/\/2, = and u- = (ho + 2/ci) 1 / 2 . We now form the ground 
state density matrix, and trace over the first ("inside") oscillator, resulting in a density 
matrix for the second ( "outside" ) oscillator alone: 

/+oo 
dxiifjo(xi,X2)iio(xi,x 2 ) 
-oo (3) 

= _ 0)1/2 exp[- 7 (x2 + 4 2 )/2 + (3x 2 x 2 ] , 

where (3 = \{u)+ — u-) 2 /(uj+ + to-) and 7 — (3 = 2u+uj-/(uj+ + u-). We would like to 
find the eigenvalues p n of p ou t(x, x'): 

/+00 
dx' pontic, x')f n (x') = p n fn(x) , (4) 
-00 
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because in terms of them the entropy is simply S = — J2 n Pn l°gPn- The solution to Eq. (4) 
is found most easily by guessing, and is 

p n = a - m n , 

(5) 

f n (x) = H n (a l ' 2 x) exp(-ax 2 /2) , 

where H n is a Hermite polynomial, a = (7 2 — /3 2 ) 1 / 2 = (cu+cj-) 1 / 2 , £ = /3/(7 + a), and n 
runs from zero to infinity. Eq. (5) imples that p ou t is equivalent to a thermal density matrix 
for a single harmonic oscillator specified by frequency a and temperature T = a/ log(l/£). 
The entropy is 

S(0 = -log(l-0-T^log£, (6) 
where £ is ultimately a function only of the ratio k\/kQ. 

We can easily expand this analysis to a system of N coupled harmonic oscillators with 
hamiltonian 

N N 

H =\^2Pi+\ Yl x i K ij x 3 ' ( 7 ) 
i=l i,j=l 

where K is a real symmetric matrix with positive eigenvalues. The normalized ground 
state wave function is 

4>o(xi,...,x N ) = 7r" Ar / 4 (detO) 1/4 exp[-x-0-x/2] , (8) 

where O is the square root of K: if K = U T K B U, where K u is diagonal and U is orthogonal, 
then O = U T K^ 2 U. We now trace over the first n ("inside") oscillators to get 

/n 
Y[ dxi ipo(xi, . . . , x n , Xn+i, . . . , x N ) 
i=l 

X Vo(^l» • • • > ^n+l> • • • ' x 'n) • ( 9 ) 

To carry out these integrals explicitly, we write 

A B 

B T C, 

where A is n x n and C is (N — n) x (iV — n). We find 

Pout(^5 x) ~ exp[— + x'-^-x')/2 + x-(3-x \ , (11) 

where x now has N — n components, (3 = |i? T A~ 1 i?, and 7 = C — /?. In general (3 and 
7 will not commute, which implies that Eq. (11) is not equivalent to a thermal density 
matrix for a system of oscillators. 
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We need not keep track of the normalization of p ut 5 since we know that its eigenvalues 
must sum to one. To find them, we note that the appropriate generalization of Eq. (4) 
implies that (det G) p ou t(Gx, Gx') has the same eigenvalues as p ut(^ 5 x ')i where G is any 
nonsingular matrix. Let 7 = V 7 D 7, where 7 D is diagonal and V is orthogonal; then let 
x = F T 7 D 1 ^ 2 y. (The eigenvalues of 7 are guaranteed to be positive, so this transformation 
is well defined.) We then have 

Pout(y,y') ^ exp[-(y-y + y'-y')/2 + y-/3'-y'] , (12) 

where j3' = 7 D 2 F/?F T 7 D 1 ^ 2 . If we now set y = Wz : where W is orthogonal and W^fl'W 
is diagonal, we get 

N 

Pout(z, z') ~ n exp [-(zf + z[ 2 )/2 + fo^ , (13) 
i=n+l 

where f3[ is an eigenvalue of . Each term in this product is identical to the p ou t of Eq. (3), 
with 7 — > 1 and (3 — * j3'-. Therefore, the entropy associated with the p ou t of Eq. (13) is just 
S = where S(£) is given by Eq. (6), and & = (3'J [l + (1 - ) 1 ' 2 ] . 

We now wish to apply this general result to a quantum field with hamiltonian 

*-l/<M« , <*>+iv**)i i l- <") 

To regulate this theory, we first introduce the partial wave components 



(15) 



where x = |x| and the Z\ m are real spherical harmonics: Ziq = Y/q, Z\ m = V2KeYi m for 
m > 0, and Z; TO = >/2Iml/ m for m < 0; the Z^ TO are orthonormal and complete. The 
operators defined in Eq. (15) are hermitian, and obey the canonical commutation relations 

Wlm{ x )^l'm'{ x ')\ = ^ll^mm>K x ~ x ') ■ (16) 
In terms of them, we can write H = J2lm ^lm^ where 



Hi m = hJ o dx \ *lm( x ) + x2 ^ ( 



9 ( <flm( x ) 



X 



2 



Z(Z + 1) 2 , s , 
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So far we have made no approximations or regularizations. 

Now, as an ultraviolet regulator, we replace the continuous radial coordinate x by 
a lattice of discrete points with spacing a; the ultraviolet cutoff M is thus a -1 . As an 
infrared regulator, we put the system in a spherical box of radius L = (A + l)a, where A 
is a large integer, and demand that <pi m (x) vanish for x > L; the infrared cutoff [i is thus 
L . All together, this yields 



tt 1 V L 2 i (a i 1\2 ( W™,3 <Plmj+l \ 2 , /(/ + !) r j 
H lm-^Z^^lm, j + {J + ^) JTTJ • —j^r, in . 



(18) 



where V?Zm,./V+l = 0> Wmj an d "'Imj are dimensionless, hermitian, and obey the canonical 
commutation relations 

[flmJ^l'm'J'] = iS ll'$mm' S jj' ■ ( 19 ) 

Thus, Hi m has the general form of Eq. (7), and for a fixed value of A we can compute 
(numerically) the entropy Si m (n, A) produced by tracing the ground state of H\ m over the 
first n sites. The ground state of if is a direct product of the ground states of each if/ m , 
and so the total entropy is found by summing over / and m: S(n,N) = J2im^lm( n ^)- 
As can be seen from Eq. (18), if/ m is actually independent of m, and therefore so is 
Si m (n,N) = Si(n,N). Summing over m just yields a factor of 21 + 1, and so we have 
S(n,N) = + l)Si(n, N). From Eq. (18) we also see that the /-dependent term 

dominates if I ^> N, and in this case we can compute Si(n, N) perturbatively. The result 
is that, for I ^> A, S/(n, A) is independent of A, and is given by 

S,(n,J\0 = £,(n)[-log&(n) + l] , (20) 

where 

n(n+l)(2n+l) 2 6 
^ (n)= 64*2(1 + 1)2 )• ^ 

Eqs. (20) and (21) demonstrate that the sum over / will converge, and also provide a useful 
check on the numerical results. 

Let us define R = (n + \)a, a radius midway between the outermost point which was 
traced over, and the innermost point which was not. The computed values of S(n, A) are 
shown for A = 60 and 1 < n < 30 as a function of R 2 in Fig. 1. As can be seen, the points 
are beautifully fit by a straight line: 



S = 0.30 M 2 R 2 , 
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where M = a . Furthermore, S(n,N) turns out to be independent of N (and hence the 
value of the infrared cutoff). Specifically, for fixed n, with n < \N ", the values of S(n,N) 
turn out to be identical (in the worst case, to within 0.5%) for N = 20, 40, and 60. The 
restriction to n < \N is necessary, since the linear behavior in Fig. 1 cannot continue all 
the way to n = N: at this point we will have traced over all the degrees of freedom, and 
must find S = 0. S must therefore start falling as R begins to approach the wall of the 
box at radius L = (N + l)a. 

Of course, similar calculations can be done for one- and two-dimensional systems as 
well. For d = 2, our introductory arguments would lead us to expect that S = nMR, 
since the relevant "area" is the circumference of the dividing circle of radius R. This is 
confirmed by the numerical results, which will be presented in detail elsewhere [5]. For 
d = 1, our arguments must break down: they would lead to the conclusion that S is 
independent of R, and this is clearly impossible. In fact, the numerical results indicate 
that S = k\ log(MR) + K2 log(pR) in one dimension; for the first time, we see a dependence 
on the infrared cutoff fx [6]. For d > 4, regularization by a radial lattice turns out to be 
insufficient; the sum over partial waves does not converge. Regularization by a full, d- 
dimensional lattice would certainly produce a finite S, but this procedure would greatly 
increase the computational complexity. 

To summarize, a straightforward counting of quantum states in a simple, well-defined 
context has produced an entropy proportional to the surface area of the inaccessible region, 
inaccessible in the sense that we ignore the information contained there. Eq. (22) is 
strikingly similar to the formula for the entropy of a black hole, S B n = ^M^A, and so may 
provide some clues as to its deeper meaning. 

I am grateful to Orlando Alvarez, Steve Giddings, David Gross, and Andy Strominger 
for helpful discussions. 

Note added: After this paper was completed, I learned of related work by Bombelli 
et al [7] . Also motivated by the black hole analogy, these authors find an equivalent result 
for the entropy of a coupled system of oscillators. They also argue that, for a quantum 
field, the entropy should be proportional to the area of the boundary; the argument they 
give is different from those presented here, and is valid only if the field has a mass m which 
is large enough to make the Compton wavelength 1/m much less than R. I thank Erik 
Matinez for bringing this paper to my attention. 
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Fig. 1. The entropy S resulting from tracing the ground state of a massless scalar field over 
the degrees of freedom inside a sphere of radius R. The points shown correspond to 
regularization by a radial lattice with N = 60 sites; the line is the best linear fit. R is 
measured in lattice units, and is defined to be n + ^, where n is the number of traced 
sites. 
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